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We show that a chiral nonUnear sigma model coupled to degenerate neutrons exhibits a ferromag- 
netic phase at high density. The magnetization is due to the axial anomaly acting on the parallel 
layers of neutral pion domain walls spontaneously formed at high density. The emergent magnetic 
field would reach the QCD scale 10^^ [G], which suggests that the quantum anomaly can be a 
microscopic origin of the magnetars (highly magnetized neutron stars). 



Magnetars and high density neutron matter. — 

The phase diagram of QCD is a mystery to be uncovered. 
Even though the problem is theoretically well-posed with 
QCD Lagrangian, so far the difficulty remains yet to re- 
veal a specific region of the phase diagram, the region at 
low temperature and high baryon number density p [1 . 

One of the promising systems where such a high den- 
sity region is realized in nature is the deep interior of 
compact stars, such as the neutron stars. Observations 
of various properties of these stars should give us cru- 
cial constraint on high density matter. In particular, 
the magnetars, which are considered to be neutron stars 
with very strong magnetic field 10^^ [G] at their sur- 
face, are of partiular interests j2| [3j. The mechanism 
for generating such a strong magnetic field is so far not 
established: Among various proposals including the the 
dynamo formation model and the fossil- field model |2], 
nuclear ferromagnetism associated with the solidification 
of the neutron star core, which was originally proposed 
right after the discovery of the pulsar [4 , is an interest- 
ing possibility to generate large intrinsic magnetic field. 
However, modern quantum many-body calculations on 
the neutron matter and asymmetric nuclear matter with 
realistic nuclear force have shown that these systems stay 
in the liquid phase at high density without having spon- 
taneous ferromagnetic transition [5 . Another interesting 
possibility is the ferromagnetism of the quark liquid in 
the central core of neutron star [6]; spin-polalized quark 
phase, similar to that in the low density electron gas, may 
be realized in a certain window of baryon density due to 
the Fock term of the gluon exchange between quarks. 

In this letter, we propose a novel mechanism which 
leads to a spontaneous magnetization of the neutron mat- 
ter, based on the non-linear chiral Lagrangian of pions 
coupled to degenerate neutrons. Two basic ingredients 
are (i) the neutron spin-density induced on a pion do- 
main wall in dense matter [7^ and (ii) the baryon number 
induced on the pion domain wall by an external mag- 
netic field through axial anomaly [8 . We show that lay- 
ers of TT^ domain walls are spontaneously generated by a 
small seed of an external magnetic field. Then, intrinsic 
magnetic field is induced from the layers which have net 



magnetization. If this mechanism takes place inside the 
core of the neutron stars above certain threshold density, 
they acquire large magnetic field and become magnetars. 

Chiral Lagrangian and pion domain walls. — We 

use the chiral Lagrangian for low energy pions and nu- 
cleons with the Weinberg parametrization [9^: 
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Here N is the nucleon field (isospin SU(2) doublet), 
is the pion (triplet), /^^ is the pion decay constant, qa is 
the axial charge of the nucleon, and is the pion mass. 
Also, = 

with D^cj) 

In [7 , neutrons are integrated with neutron chemical 
potential, and tt^ domain wall solutions were studied. 
Here we generalize it to include the charged pion and 
obtain the in-medium chiral Lagrangian up to O(p^); 
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with i = 1, 2, 3. We have defined the charged pion = 

^i+z^2, and D^cj)^ = + Z(5^oMi^+)/(l + |0| V^/^), 

where /ii is the isospin chemical potential defined as the 
difference /ii = /ip — /in- The correction due to the 
background neutrons is in the coefficients a, /3, a, P (and 
71,2,3,4), which are equal to the unity (zero) in the ab- 
sence of the nucleons. They are given by one- loop calcu- 
lations (see [7 for a and /3) 
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FIG. 1: Plots of P (rigid lines) and /S (dashed lines) as a func- 
tion of /in. P decreases monotonically as the density grows. 
Here we used 940 [MeV], 93[MeV] and qa = I 

(with which /3 vanishes at /in I.OSttin corresponding to 
Pn 0.23 [fm~^]). In general in dense matter and /^r may 
vary, so the plot should be understood only qualitatively, a 
and a do not vary much in density. 
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Here x = jin/'^^i and we consider the pure neutron mat- 
ter for simplicity, so that /ip = . The neutron density 
is related to the Fermi momentum as pn = /cp/(37r^) with 
/cf = fJ^n ~ "^^^ correction is taken care of at 
its leading order. An important feature of the correc- 
tions is that /5 is a monotonically decreasing function of 
the neutron density, and vanishes at a certain density 
(Fig.[l]). 

A classical solution of Eq.Q is a domain wall of the 
in-medium neutral pion [7], 
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Note that this tt^ domain wall interpolates the vacua = 
and = 27r, where tan | = 1 01 72/7^. Interestingly, the 
domain wall can reduce its weight in the neutron matter: 
70 stays positive at all x, while /3 approaches zero as 
shown in Fig.l, so that the tension £ / S of the domain 
wall is significantly reduced for ^5 ^ 0, 



£/S = 8v^/>., 
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where S is the domain wall area. We will show that the 
parallel layers of the domain walls would populate at high 




FIG. 2: The multi domain wall solution written by 0{x^). 
Each wall interpolates adjacent vacua, = 0, 27r, 47r, • • • . For 
the exact solution, see for example ^lOj . 



baryon density due to the reduction of the tension and 
the accumulation of baryon number on the domain wall. 

Emergent magnetic field on pionic walls from ax- 
ial anomaly. — Let us first discuss how the spon- 
taneous magnetization occurs in high density neutron 
matter. There are three steps for this to happen: (i) 
The in- medium domain wall becomes light, and at the 
same time it acquires finite baryon density due to ax- 
ial anomaly under an external magnetic field. Then the 
system with a domain wall becomes energetically favor- 
able than the uniform neutron matter above a certain 
density, (ii) The domain wall is magnetized due to the 
spin alignment of surrounding neutrons, so that it cre- 
ates spontaneous magnetization which enlarges the orig- 
inal magnetic field, (iii) The enhanced magnetic field 
creates more domain walls. The cycle (i)^(ii)^(iii)^(i) 
is repeated and leads to a stable configuration with many 
parallel layers of thin domain walls with a high magnetic 
field. All the spins are aligned, so the system is ferro- 
magnetic. 

We now explain the mechanism of inducing the baryon 
charge following Ref.[8 . In a background constant mag- 
netic field, the QCD axial anomaly term in the chiral 
Lagrangian reads 
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where Aq is the temporal component of a gauge po- 
tential for the baryon number symmetry, Bs is the back- 
ground magnetic field along x^, and U = cos 6>+iT-0 sin 
with (f) = cj)/\cf)\. Since the pion-dependent part can be 
evaluated (for = 02 = 0) as 
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we immediately see that the domain wall, which interpo- 
lates = and ^ = 27r, can obtain a baryon charge per 
a unit area 1 8 , 



Nb/S = 6^3/2^. 
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If the domain wall is not parallel to the magnetic field, 
the induced baryon charge is reduced to Bihi where n is 
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the unit vector perpendicular to the domain wah. Note 
that the formula Eq.([7|) is valid even for /3 / 1. 

Combining Eq. ^ with Eq. Q , the domain wall energy 
per a unit baryon charge is 



(8) 



The system with domain wall is more favorable than the 
uniform neutron matter when the domain wall energy 
per baryon becomes smaller than the neutron chemical 
potential [8 , i.e. Bs > V f^lo x (167r/^m7^//in) ~ V f^lo x 
10^^ [G]. Here, the factor /3, which was not taken into 
account in [8 , is important. Since f3{p) is a monotonically 
decreasing function of pn, an adiabatic increase of the 
density inevitably hits the critical value of /3 at which 
the domain wall is created. 

Once the pionic wall is formed, neutron spins on the 
wall align in the direction perpendicular to the domain 
wall [7] . The spin density of the neutrons is a spatial part 
of the axial current j^^"^ = (V^nTiTsV^n)- It was evaluated 
in 0, in the same approximation, as 



S3/5 = 27r(/3-l)/2, 
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which is the third component of the neutron spin density 
per a unit area of the domain wall. 

We note here that neutrons have a magnetic moment 
/i = gess/^m^ where g ~ —3.8 is the neutron ^-factor. 
So the total magnetic moment (which is the magnetiza- 
tion M) per a unit volume at /3 ~ is 
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where d is the separation between adjacent domain walls. 
Therefore the domain wall phase is ferromagnetic. 

The magnetization M is larger for smaller separation 
d among the domain walls. This d is intimately related 
to the induced baryon density due to the domain wall, 
and in fact this is a driving force for developing strong 
magnetic field. From Eq.Q, we know that the averaged 
baryon number density induced by the domain walls is 
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The total baryon number density p is given as a sum 
of pdw and the remaining neutron density pn- (This pn 
cannot vanish, since the (3 correction needs background 
neutrons.) Combining Eq.(pT]) with Eq.(lO), we obtain 
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Once the background Bs becomes larger, the domain wall 
energy cost ([8| becomes smaller. So the domain wall 
is created easier, and the domain wall separation d be- 
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FIG. 3: The topological path of the neutral pion domain wall. 
The sphere represents the target space of the sigma model. 
On the S^, the thick line with an arrow represents the domain 
wall solution (|3|. It rounds the sphere, but topologically triv- 
ial. The plane beneath the sphere is the parameterization 
space of 03 and 01,2. There is a one-to-one correspondence 
between a point on the sphere and a point on the plane. 



increases and helps the Bs to increase. So, this system 
has a self-enhancement mechanism of the magnetic field. 
The equilibrium can be reached at ^3 = M which is our 
critical induced magnetic field. 



comes smaller. Then from Eq.(lO), the magnetization M 



B3 = V27r2|5|/2pdw/m^ - 3 x lO^^ [G]. (13) 

supposing a typical value for p ~ p^w ^ Pn- The value 
of the magnetic field is quite large {^/Bs ~ 10^ [MeV]) 
and close to the QCD deconfinement scale around which 
our approximation breaks down [15 . The magnetization 
is expected to stop increasing somewhere before reaching 
this value. 

In addition to the neutron spin alignment, the WZW 
term ([5| itself may provide a magnetization [13] of the 
same sign. Details will be in our forthcoming paper [M]. 

Finally we briefiy comment on the stability of the do- 
main walls. Our domain wall is topologically trivial, be- 
cause the target space of the non-linear sigma model is 
SU (2) on which the two vacua ^ = 0, 27r are the 

same point (south pole): See Fig. |2] and Fig. [sj There- 
fore, the domain walls can be created spontaneously, but 
on the other hand, they can decay through the fiuctua- 
tions along 0i^2 directions. Indeed, it was shown that the 
domain wall is stable only for Bs > 10^^ [G] at /3 = /3 = 1 
by analyzing its local and global stability [8 . As we have 
shown after Eq.(IT3|, the above condition is replaced by 
Bs > y/Wlo X 10^[G]. This implies that the global sta- 
bility is guaranteed for smaller Bs when /3 ^ 0. The 
analysis of the local stability is more involved especially 
where /3 > and < (see Figjl]): In this case, the 

charged pion condensation 0+ = aexp (^—ipT^t-\-ik'X^ 

should be considered together with the tt^ domain wall. 
We leave the analysis to our future work. p!6l 

Implication to magnet ars. — A schematic picture of 
the core region of the neutron star is shown in Fig. [4] 
Only at the core region, because of the high density (or 
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Layers of neutral 
pion domain walls 



FIG. 4: A schematic figure of the neutron star with domain 
wall layers at the core. Scales should not be taken seriously. 



rather to say the large value of the neutron chemical po- 
tential), the domain walls are present. At the boundary 
of the domain walls, neutrons drip from the wall bound- 
ary so that the total baryon number is conserved [TT]. 
As the domain wall layers, the ferromagnetic region, are 
present only at the core of the neutron star, at the sur- 
face of the neutron star the magnetic field does not reach 
the value in Eq.(13), but it would be enough strong to 



explain the magnet ars. 

The magnetic field at the surface of the neutron star is 
smaller than that of the domain wall core, as ^surface = 
^core(^core/^Ns)^ at the North pole, i^core IS the radius 
of the domain wall core (assumed to be spherical and 
to have a homogeneous ferromagnetism inside) , and R^s 
is the radius of the neutron star. The average of the 
magnetic field magnitude on the neutron star surface is 

5; 



average 



0.695surface- The corc-radius dependence of 
the surface magnetic field is shown in Fig. |5] If the core is 
sufficiently small such as i^core/^NS ^ the surface 

magnetic field may reduce to ^surface ~ O{10^^) [G]. 

The mechanism suggests that there are two kinds of 
neutron stars: one which reaches the critical density and 
has the domain wall layer structure, and the other which 
does not have it. The former would have a strong mag- 
netic field but the latter would not have it. It is inter- 
esting that the recent data [3 of the magnitude of the 
magnetic field on the surface of the neutron stars show 
two categories, magnetars and the others. 

It is important to construct more realistic models with 
nuclear forces, as our model uses free neutrons. For ex- 
ample, the neutron superfluidity can co-exist with the 
domain wall, since at higher densities the spin-aligned 
neutron pairing is known to be favored. Furthermore, 
structure of the solitonic core of neutron stars would in- 
fluence the equation of state, and may be sensitive to 
the mass-radius map of the neutron stars. If the orien- 
tation of the solitonic core is different from the rotation 
axis of the neutron star, it would be a source of gravita- 
tional waves. All details need to be explored, to match 
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FIG. 5: The magnetic field at the neutron star surface. i?core 
is the radius of the domain wall layer core (assumed to be 
spherical), and i?surfce is the radius of the neutron star. At 
the core, the critical magnetic field (13) is assumed. For the 
I.4IM0 neutron star with the standard APR equation of state 
[12^, Rcove/Risis =0.1 corresponds to the critical density of the 
pionic wall formation 3.2po- 



the observations of the neutron stars. We hope that our 
mechanism may survive various corrections, and explain 
observations. 
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